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It is shown, how the normalization sum rules of the spin-parallel and spin-antiparallel pair 
densities of the homogeneous electron gas become sum rules for the scattering phase shifts of 
the Overhauser two-body wave functions (geminals), with which the pair densities have been 
successfully parametrized recently. These new sum rules relate two-body quantities to a one-body 
quantity, namely the asymptotics of the Overhauser geminals to the momentum distribution. 
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The (spin-unpolarized) homogeneous electron gas 
(HEG) is an important and widely used model for 
the phenomenon called electron correlation, cf. e.g.El. 
Some of its details are hidden in the reduced densi- 
ties. For the HEG ground state such characteristic 
quantum-kinematic quantities are (i) the momentum dis- 
tribution n(k) as the typical 1-body quantity (recently 
parametrized in terms of a special convex functional) and 
(ii) the pair density (PD) g(r) as the_typical 2-body 
quantity (recently parametrized in Refs.B-El). With spin- 
resolution one has to distinguish the PDs g^(r) and 
g-\l(r) for electron pairs with an interelectron distance 
r = |tj — r 2 \ and with parallel and antiparallel spins, 
respectively. They describe the probabilities of finding 
in the neighborhood of an electron (the 'blue' electron) 
another electron with the same respectively the opposite 
spin. So it holds g^{r) > and g-\i(r) > 0. For large r it 
is <7i"|-(oo) = 1 and g^i(oo) = 1. For r = ('on-top') it is 
fl , tt(0) = due to the Pauli repulsion and <7f;(0) < 1 due 
to the Coulomb repulsion. These short-range deviations 
of g-\-\(r) and gn(r) from their asymptotical values 1 are 
called Fermi hole and Coulomb hole, respectively. In be- 
tween they-, show shell structures. They are normalized 
as (cf. e.g.El) 



p / d 3 r[l-g n (r)}=0 (1) 



d 3 r[l- m (r)] = 2, 



with p = electron density. The spin-summed PD g(r) = 
^ [ffTT ( r ) + 9ll( r )] i s therefore normalized as p J d 3 r [1 — 
g(r)] — 1, what is called perfect screening sum rule or 
charge neutrality condition. Here and in the following 
wave lengths and momenta k are measured in units of the 
Fermi wave length fcp = a = and lengths 

(like r or p^ 1 ^ 3 ) in units of fcp 1 . Therefore p = -^j 
and pd 3 r — a 3 d(r 3 ). The momentum distribution and 
the PDs depend parametrically on r s , the radius of a 
sphere containing in the average one electron, measuring 
simultaneously the interaction strength with. r s = cor- 
responding to the ideal Fermi gas. In RefsJj-13 the PDs as 
functions of the separation r and of the parameter r s have 



been fitted to the data of quantum Monte-Carlo calcu- 
lations. With the help of (such available) PDs, particle- 
number fluctuations in fragments of the system can be 
discussed.0 The partitioning into a generalized Hartree- 
Fock term and a connected remainder 



9n(r) = 9n(r) h n (r), gff(r) = 1 - |/(r)| a , 



HF/ 



9n(r) = 9W(r)-h n (r), gff (r) = 1 



HF/ 



(2) 



defines the cumulant PDs h^(r) and /ifj(r). f(r) = 
{1/N) J^fc n(k)e— is the dimensionless 1-body reduced 
density matrix. /(0) = 1 expresses the normalization of 
n(k). The hole normalizations (1) imply 



<Prh rt (r) = 2c, P d 3 rh n (r)=0, (3) 



where c —(2/N) J^k n (k)[l — n(k)] is christened Lowdin 
parameters (measuring the correlation induced nonidem- 
potency of n(k), thus vanishing for n(k) — > — 
0(1 — k)). It is easy to show that 

.g T H T F (r)=2 5 _(r), g™(r) = g + (r) + g_(r), (4) 

define a singlet term g+ F (r) and a triplet term g^fr) 
with 



<?± F M= V <jf(kr)>, L=(l, mi ), 



(5) 



where ± refer to even and odd I, respectively, and 
— J2i(2l + 1) • • ■ • The fc-average is defined 
as < ••• >= (2/AOEfcM(*0---, (VN)J2Mk) = 1, 
p(k) = (2/N)J2Kn(k 1 )n(k 2 ), where k ia = \\K±k\. 

Note (2/N) Efe • • • = XT rf ( fc3 ) ' ' ' • Thus ' 3fc2 M( fc ) is the 
probability of finding two electron momenta k x and fc 2 
with the half difference k. Because of 



d 3 r 
d 3 r 



9T(r) \ 
9±{r) - ^ 



= ±(l-c), 
= ±1, 



(6) 
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the corresponding cumulant PDs h±(r) = g± — g±( 



are normalized as 



p J d 3 rh±(r) = =Fc. 



(7) 



g±(r) and g± F (r) approach 1/2 for r — > oo, thus 
fe±(oo) = 0. Note fl (r) = [g+(r) + 3<?_(r)]/2. 

No correlation (corresponding to r s = 0) means 
ra(Jfc) -» n<°)(fc) = 9(1 - Jfe), thusrtt*^) ^ (0) O) = 
4(1 - fc) 2 (2 + fc)0(l - fc), cf. RefJla, and h±(r) = 0. 
In this ideal case the identities 



sin 2fcr 
< 2kr > 



3ji(r) 



sin 2fcr 
2kr 

2 



ensure the hole normalizations (1) to be obeyed. 



(8) 



Correlation (corresponding to r s ^ 0) means n(k) ^ 
6(1 - k) and h(r) ^ 0. Recently the HEG-PDs 
for this case have been successfully parametrized in 
terms o£ Qyerhauser geminals (= 2-body wave functions) 
Ri(r, fc)B-li3, namely as 



9±{r) 



< Rf(r,k) > 



(9) 



So the normalization (7) can be written as 



P J d 3 r^ L < [Rf(r,k) -if(fcr)] >= ±c (10) 

The question arises whether the lhs can be expressed 
in terms of quantities characterizing the differences be- 
tween the Overhauser geminals Ri (r, k) (for interacting 
electrons) and the Bessel functions ji(kr) (= Overhauser 
geminals for noninteracting electrons). Here it is shown 
that these characteristic quantities are the phase shifts 
r]i(k), which describe the large— r asymptotics of the 
Overhauser geminals Ri (r, k) according to 

RAr 1 k) — > — sin ( kr — I — h r/i(k)) for r — > oo. 
kr V 2 / 

Namely, within the Overhauser approach a pair of elec- 
trons with momenta k x and k 2 moving in the HEG is 
considered. Their center-of-mass motion is decribed by 
exx>(i KR ) , where K = fc x + k 2 * s the total momentum and 
M. = fci + r 2 ) / 2 is the center-of mass coordinate. Their 
relative motion is described by Rl(l, &) = Ri(r, k)YL(e r ) 



with the relative coordinate r 



r, and the half 



momenta difference k = 1 1 fcj — fc 2 1 ■ So the total 2-body 
wave function describing the pair is exp(i if R )Rr,(r, k). 
The Ri (r, k) are the solutions of the radial Schrodinger 
equation for the relative motion 



r dr 2 



1(1 + 1) 



Ri{r,k) = (12) 



with vf s (r) 



as an effective interaction po- 



tential. The dots indicate the potential of the physically 
plausible screening cloud around each electron (result- 
ing from exchange and correlation) , possibly different for 
"+" (even I) and "-" (odd I). Note that Eq. (|l|) is 
the remainder of an effective 2-body Schrodinger equa- 
tion and that Ri (r, k) is (part of) a 2-body wave func- 
tion, a so-called geminal. Because the effective inter- 
action potential v^ s (r) is repulsive, the solutions of Eq. 
(|l2|) are scattering states, characterized by the afore men- 
tioned scattering phase shifts r/i(k). Therefore the inte- 
grals J d 3 r jf(kr) and / d 3 r R 2 (r, k) do pot exist, but (the 
following proof is analog to that in Ref.Eil) understanding 
the lhs's of Eq. (9) as / d 3 r • • • = liniB-^oo f r<R d 3 r ■ ■ ■ 
and using the identity 



d 1 (dR*r dR L 



d 2 R, 



+ c.c. 



R* l Rl (13) 



dr 2 V dr d(k 2 ) * L drd(k 2 ) 
(this follows from Eq. (p~2h) and the Gauss theorem 



r<R 



d 3 r-^[e r f(r)} 



ArrR 2 f(R), 



it turns out 



d 3 r[R 2 (r,k)-jf(kr)] 



(14) 



(15) 



r<R 



^R 2 fdR l dR l &Ri 31 djt | d 2 ji \ 
2k \ dr dk drdk dr dk drdk ) 



r=R 



The rhs gives with Eq. (|ll| ) the sum Ai(k) + Bi(k), where 

(16) 



Bi(k) 



k 3 



[sin (2kR - In + 2r)i(k)) - sin(2fcii - In)] 



These expressions have to be /c-averaged (with p(k)) and 
L-summed. p(k) has the properties p(0) — 8(1 — c), 
n(oo) = 0, p'(k) < 0. The discontinuity of n(k) at k = 1 
(quasiparticle weight Zp) makes p"(k) also discontinous 
at k = 1 and the correlation tail n(k > 1) ^ causes a 
corresponding correlation tail p(k > 1) ^ 0. The phase 
shifts are gauged as ??z(oo) = 0. Because they originate 
from a repulsive potential, it is t]i(0) = 0. (For an attrac- 
tive potential the Levinson theorem would say rji (0) = nn 
with n =^-pumber of bound states in the I— channel, cf. 
e.g. Ref.EZl orta, p. 133.) In addition to this, for repulsive 
potentials vf s (r) with finite range it holds rji(k) ~ k 2l+1 . 
The fc-average < • ■ • > of Ai(k) yields after partial inte- 
gration 



<A ; (fc)>-6vr/ dk[~p'(k)] m (k). 



(17) 
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The other term gives < Bi(k) >i + < Bi(k) > 2 with 

/>oo 

<B ; (fc)>i = / dfc&i(fc)sin(2fci?-/7r), 
Jo 
3n 

h(k) = —n(k)[cos2r n (k)-l], 

poo 

<Bi(k)> 2 = - dkb 2 (k)cos(2kR~ In), 
Jo 



3n 

h(k) = — li{k) sin2?7;(fc). 



(18) 



Successive partial integrations yield (semi-convergent) se- 
ries starting with 

< B,{k) >x= - 6?(1+)] cos(2i? - Ztt) 

+ ^3^(0)COB/7T + .-- , 



the natural geminals, which parametrize/diagonalize 
the 2-body reduced density matrix, are presumably 
the solutions of a 2-body Schrodinger equation with a 
non-local effective interaction potential, similarly as the 
natural orbitals, which diagonalize the 1-body reduced 
density matrix, follow from an 1-body Schrodinger 
equation with a non-local effective 1-body potential. 

Whereas Eq. (|2^) is valid for a 2-body quantity of 
a uniform system (namely the HEG-PDs parametrized 
by scattering-state Overhauser geminals), the well-known 
Friedel sum rule for point defects in metals refer to a 
1-body quantity of a non-uniform system (namely the 
charge distribution or screening cloud around the impu- 
rity). This cloud perfectly screens the impurity (in a 
spatially oscillatory manner - Friedel oscillations). The 
resulting impurity potential generates scattering states 
with phase shifts T}i(k), such that the 1-body sum rule 



<Bi{k) > 2 



(2R) 



;b' 2 (0) cos In 



(19) 



which vanish for R 
rules 



oo. So the 2-body phase-shift sum 



dk [-//(*)]»»(*) = ±c 



(20) 



follow. These new sum rules are conditions for the ef- 
fective interaction potential i> efr (r), which generates the 
Overhauser geminals Ri(r,k). Whether-the phase shifts 
of the available Overhauser geminalsB-E3 really obey 
the sum rules (20) has to be checked. For r s = (no 
correlation) the Lowdin parameter c and the Overhauser 
phase shifts rji(k) vanish. 

The new sum rules ( p0| ) are relations between the 
one-body momentum distribution n(k) encoded in /J.(k) 
and c and the asymptotic behavior of the PD for large r 
as described by the phase shifts r/i(k) of the Overhauser 
two-body wave functions Ri (r, k) . This may be consid- 
ered as complementary to the large-A: asymptotics of 
n(k), which is determined by the on . top .PD according 
to n(k -> oo) ~ Su(°)/ fcS + O(l/k 10 )HM A qualitative 
discussion of the mutual relation between n(k) and g(r) 
based on the virial theorem is in RefBa 



(21) 



holds. The lhs is the excess charge or the valence 
difference between the impurity and the host and the rhs 
is the total number of electronic states induced by the 
impurity. For the details of this subject cf. e.g. Refs.lia, 
p. 167, 469 oilij, p. 44. For the generalisation of the 
Friedel sum rule to non-spherically symmetric scatterers 
cfEj. For the Sugiyama-Langreth neutrality rSijm rule 
of half-space jellium and its generalisation cf.Eil andE3, 
respectively. 

Note the essential difference between Eq. (|2~i|), where 
an external charge causes a perturbation, and Eq. (|2C|), 
where the 'blue' electron, which distorts its surrounding, 
is an internal charge belonging to the quantum many- 
body system with the consequences of (i) exchange 
(singlet/triplet or ± symmetry) and (ii) of the k average 
<•••>. 
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